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THCIINICAL MiMORANDUM 


A STUDY INTO THE LOSS OF LOCK OF THE SPACE TELESCOPE 
FINE GUIDANCE SENSOR 

SUMMARY 


This report documents tlie results of a study into the loss of lock phenomenon associated with 
the Space Telescope (ST) Fine Guidance Sensor (FGS). The primary cause of loss of lock has been 
found to be a combination of cosmic ray spikes and photon noise due to a 14.S Mv star. The probability 
of maintaining lock versus time is estimated both for the baseline FGS design and with parameter changes 
in the I’GS firmware which will improve the probability of maintaining lock. The parameters varied are 
changeable in-flight from the ground and hence do not Impact the design of the FGS hardware. 


I. INTRODUCTION 


The Space Telescope (ST) is an unmanned, Farth orbititig astronomical observatory to be inserted 
into a oOO km circular orbit by the Space Shuttle in the mid-1980’s, The ST will have the capability of 
pointing at celestial targets with an accuracy of 0.01 arc-sec and a stability of 0.007 arc-.sec RMS for 
observation times up to 10 hours [1 ]. This is accomplished by the onboard pointing control s> item 
(PCS) winch includes a system of rate gyros, reaction wheels, an onboard digital computer, and a fine 
guidance sensor (FGS) [21. The fine guidance sensor is basically a system of three star trackers located 
at the focal jilane of the telescope. Two out of the three star trackers are normally employed to generate 
vehicle attitude information whicli is used to periodically update that derived from the rate gyro outputs. 
Fach star tracker acquires and tracks known guide stars as dim as 14,5 Mv. Tracking is accomplished by 
means of a fine lock loop, internal to each star tracker, which consists of a two u.\is Koesters prism 
interferometer, four photomultiplier tubes, digital control logic in the fine guidance electronics (FGF), 
and two star selector servos (glmbai rate servos), Optical elements are mounted on the star selector 
servo (SSS) shafts which deflect the guide star image in such a way so as to compensate for vehicle line 
of sight motions in order to keep the guide star at null in the field of view of the interferometer. Out- 
puts from the photomultiplier tubes and outputs from the gimbal encoders mounted on the SSS shafts 
are routed to the PCS digital computer every 0.025 sec and used to update the computed veliide attitude 
derived by integrating the rate gyro outputs. These updates take place once per second. However, 
disturbances to the FGS fine lock loop caiv cause the FGS to lose lock on the guide star it is tracking. 

To date, the principal cause of loss of lock has been found to be a combination of cosmic ray spikes 
striking the PMT photocathode and photon noise due to a 14.5 Mv star. 


This report presents the results of a studj' 


LJ f/l * 


a) Analyze the loss of lock phenomenon associated with the FGS fine lock loop. 

b) Predict the probability of maintaining lock with the FGS, calculated based on loss of lock 
results from a single a.xis digital simulation of the ST PCS/FGS. 

c) Offer potential solutions for improving the probability of maintaining lock wliicli do not 
require redesign of the FGS. 


II. DESCRIPTION OF THE FGS FINE LOCK LOOP 


original Rage is 

OF POOR QUAUr/ 


The primary purpose of the FGS fine lock loop is to keep the guide star image at null in the field 
of view of the two axis Koesters prism interferometer. A simplified functional block diagram of one 
axis of the fine lock loop is shown in Figure 1. The outputs of the interferometer are the photons from 
the guide star which, in turn, strike the photocathodes of the PMTs, two pet axis. The photons generate 
photoelectron counts in the PMTs and these are integrated over 0,025 sec intervals and read out by the 
FGE. These are inputs into the calculation of Q, the normalized fine error signal. The transfer charac- 
teristic for the interferometer from the input to the interferometer, denoted by to the normalized 

error signal, Q, is shown in Figure 2, both for the ideal theoretical curve and the simplified model used 
in the simulation. Because of the uncertainty that exists as to the location of the tails of the actual 
manufactured interferometers, the simplified model used in the simulation goes to zero at ±0.047 arc-sec. 

Returning to Figure 1, the output of the Q computation is limited to ±1 and this result is input 
into the linearization algorithm which; 

a) Transforms the normalized error signal into object space. 

b) Helps to linearize the interferometer transfer characteristic in the region between ±0,01 arc- 
sec for Pj. 

Nominal values for K^, K^, and K 3 in the linearization algorithm are given in Figure 1, The output R of 
the linearization algorithm is then transformed mto gimbal space by the optical lever arm Kqq after a 

sign inversion, and then tWs result is input into the fine lock loop proportional-integral (P-I) control law 
which generates gimbal rate commands to the star selector servo (SSS) every 0,025 sec. The gimbal rate 
commands to the SSS result in a gimbal rate and hence a change in shaft angle which, in turn, moves 
the star image by virtue of the optical elements on the SSS shaft in the optical train on the FGS. The 
action of the fine lock loop compensates for line of sight motions, such as those generated by rate gyro 
noise and reaction wheel induced oscillations, to try to keep the star image at mterferometer null and 
hence in the field of view of the interferometer. However, because of the finite bandwidth of the fine 
lock loop (gain crossover frequency is about 4 Hz), it is possible for disturbances to move the star image 
outside the field of view of the interferometer. If there are no immediate subsequent disturbances to 
the fme lock loop which drive the star image back into the interferometer field of view, then loss of 
lock occurs and telescope pointing is affected. 


ill. AN ANALYSIS OF THE LOSS OF LOCK PHENOMENON 


As stated previously, loss of lock of the FGS on a guide star it is tracking occurs when disturb- 
ances to the fme lock loop cause the star image to move outside the field of view of the interferometer 
(±0.047 arc-sec in the simulation model) and there are no immedi. subsequent disturbances which 
drive the star image back into the interferometer field of view. Four disturbances to the fine lock loop 
have been simulated to date: 

a) Line of sight motion due to rate gyro noise 

b) Line of sight motion due to reaction wheel induced oscillations 
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Figure 1. FGS fine lock loop functional block diagram. 














Figure 2. Interferometer transfer characteristic. 
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c) Photon noise 


d) Cosmic ray spikes. 

It was found that rate gyro noise has little effect on loss of lock. Simulation runs were made 
with line of sight motion due to reaction wiieel induced oscillations at a frequency of 50 Hz (tliese dis- 
turbances have been predicted to be predominantly in the iTcquency range of 50 Hz to 60 Hz) and a 
zcro-to-peak ajnplitude equal to 0.01 arc-sec (it is predicted that the requirement of 0.003 arc-sec RMS 
due to reaction wheel induced oscillations will be satisfied), and it was found the oscillations actually 
alleviated the loss of lock problem because their effect was to extend th*' field of view of the inter- 
ferometer. li turns out that loss of lock is caused primarily by the combination of cosmic ray spikes 
and photon noise due to a 14,5 Mv star. At a 13 Mv star, loss of lock does not occur with or without 
cosmic ray spikes. Also, the FGS has no requirement to track stars dimmer tlian 14,5 MV for PCS 
control. 


\ 


To understand what combinations of cosmic ray spikes and photon noise from a 14.5 Mv star 
cause loss of lock, a simplified analysis was performed with a .single axis simulation of the fine lock loop, 
Initially, all line of sight motion, as well as photon noise, was set equal to zero so that the only source of 
disturbance to the fine lock loop was a single cosmic ray spike, affecting one 0.025 sec sample from the 
output of PMT A. The amplitude of the cosmic ray spike was such that it produced 300 photoelectron 
counts in PMT A as shown in Figure 3. Also shown in Figure 3 is the response of the fine lock loop, 
in terms of versus time, to this disturbadcs.. It can be seen from tlic fine lock loop response that: 


a) 0*025 arc-sec, wliich is well within the field of view of the interferometer (i 0.047 

arc-sec) and lienee loss of lock did not occur 


b) liS'i’Ijiiax 0 ‘^curred between 0.05 sec and 0,1 sec 

c) The response settled out well within 1 sec of time. 

Next the amplitude of the cosmic ray spike was varied and recorded. This was done 

for = 0.625 (nominal) and K 3 = 0, The results were plotted and are shown in Figure 4, Several con- 
clusions can be drawn from Figure 4: 

a) The value of li^xlj^ax ii'icreases with increasing amplitude for the cosmic ray spike until the 

amplitude is in the neighborhood of 60 pliotoelectron counts. This is because with a 14,5 Mv star, Q 
becomes 1 when the amplitude of the cosmic spike is about 60 counts and the Q limiter saturates at 
±1. Hence, tlie output of the Q limiter is ±1 for amplitudes of cosmic spikes greater than about 60 
pliotoelectron counts. 

b) Consequently, IjSxljuax O'-h' ‘"md at a value of 0.025 arc-sec for K 3 = 0.625 and 0.014 
arc-sec for K 3 = 0. Since the field of view of the interferometer is ±0.047 arc-sec, one cosmic spike (or 

one photon noise spike for that niattr r) will not cause loss of lock. In fact, one would predict that two 
spikes in the same tube large enough to saturate the Q limiter would be needed for loss of lock to occur 
with IC 3 = 0.625, while four spikes would be needed with K 3 = 0, Hence, one would expect loss of lock 

to be improved by reducing K 3 belov its nominal value of 0.625, 
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Figure 4. amplitude of cosmic ray spike. 





To test the theory that two cosmic ray spikes could cause loss of lock, the response of the fine 
lock loop to two cosmic ray spikes was simulated. The results are shown in Figure 5. The cosmic spikes 
were large enough to saturate the Q limiter and had the same polarity (were inserted in the same PMT), 
When they occurred back to back, or else one 0.025 sec sample apart, loss of lock occurred. When they 
were two 0,025 sec samples (or O.Oi. sec) apart, loss of lock did not occur. It can also be seen in the loss 
of lock cases that once loss of lock occurs, the star image is driven away from null at a constant rate. 

This is due to the fact that the integral path in the fine lock loop control law sums the two cosmic spF 
so that its output is a constant resulting in a constant rate command driving the star hnage away from 
null. One would suspect that reducing the value of the gain Kj 4 in the integral path would improve the 
ioss of lock problem, also. 

It should be pointed out that it sliould be possible to find combinations of three or more smaller 
amplitude noise spikes which will cause loss of lock, This was not investigated, however. What has been 
established is that one noise spike (cosmic ray or photon) will not cause loss of lock regardless of ampli- 
tude because of the Q limiter in the FGE, Two noise spikes of the same polarity (in the same tube) 
large enoug<’< to -saturate the Q limiter and close together will cause loss of lock. 

To further understand the nature of the loss of lock problem, simulation runs were made with a 
single axis simulation of the PCS/FGS under more realistic conditions with cosmic spikes of various 
amplitudes. In these simulation runs, the PCS/FGS was configured in what is known as Configuration 2 
wherein the FGS outputs (PMT counts and gimbal encoder readouts) every 0.025 sec are used to update 
once per second the vehicle attitude computed onboard from the rate gyro outputs. This is the normal 
rnede for fine pointing. It was assumed that one rate gyro had failed and of the remaining five, the best 
four were chosen so as to minimize rate gyro noise. There was no reaction wheel induced jitter. The 
FGS was tracking a 14.5 Mv star and the two axis photon noise equivalent angle was 0.0046 arc-sec. 
Cosmic spikes were inserted in PMT A every 2 sec. 2 sec was chosen based on the 1 sec settling time of 
the fine lock loop. Hence, the effects of one cosmic spike settled out before the next arrived, which is 
what happens in tJie actual situation as shall be seen later. Finally, the simulation run times were 5400 
sec long, The amplitudes of the cosmic spikes in a given run were held constant, but varied from run to 
run. It was ohseiwed whether loss of lock occurred in 5400 sec. The results arc shown in Table 1. 

The results indicate that the FGS maintains lock until the amplitude of the cosmic spikes are 52 photo- 
electron counts or more. Hence, it is seen that loss of lock does not occur unless cosmic spikes are pre- 
sent and their amplitude exceeds approximately 50 counts. Should cosmic spikes greater than 50 counts 
exist in the actual situation, one would intuitively expect the probability of maintaining lock to be 
related to the number of these events that occur. 


In a separate study, summarized in Table 2, it was found that the mean time between cosmic 
spikes greater than 50 photoclectron counts, hitting any one of eight PMT’s (four PMT’s per star tracker 
with two normally used for control) was 673 sec with a pulse amplitude discriminator (PAD) setting of 
80 photoclectron counts in the FGS and 1293 sec with a PAD setting of 40 photoelectron counts, ^ 

The PAD is a setthig in the electronics at the output of each PMT wWch discriminates against some of 
the cosmic spikes or photon noise spikes which exceed the setting. Unfortunately, the PAD does not 
discriminate against all the cosmic spikes which exceed the threshold setting and its ability to do so 
depends on the direction of the cosmic rays relative to the PMT. It has been predicted that 40 is a 
consei-vative value for the PAD setting at the beginning of life and 80 at the end of life, the change 
caused by degradations in the electronics with time. These results will be used in Section IV of this 
study to estimate the probability of maintaming lock with the FGS. 


1 . This is the result of an in-house study by William B. Chubb of NASA-Marshall Space Flight Center 
in February, 1983, 
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IV. THE PROBABILITY OF MAINTAINING LOCK WITH THE BASELINE SYSTEM 


It has been established that loss of lock occurs when the fine lock loop tracks a 14.5 Mv star 
and cosmic ray spikes exceeding 50 photoelectron counts are present. However, time to loss of lock, 
as it turns out, is not deterministic, but depends on the probability of getting a large cosmic ray spike 
surrounded by the right photon noise profile, To determine the probability of maintaining lock as a 
function of time, the preferred approach would be to simulate the cosmic spikes at the frequency they 
are expected to occur in the actual situation, wliich is once every 673 sec or 1293 sec, depending on the 
PAD setting. The simulation would be run until loss of lock and then this procedure would be repeated 
a number of times varying the photon noise profile from run to run to get a random sample of the time 
to loss of lock. From tliis, the probability of maintaining lock versus time could be estimated. However, 
simulation run times would be intolerably long to do this. Consequently, the scheme employed was: 

a) Insert cosmic spikes every 2 sec in PMT A in the simulation. (2 sec was chosen to be as small 
as possible to save computer run time but larger than the 1 sec settling time of the fmc lock loop to a 
single cosmic ray spike. Hence, the effect of one spike will have settled out before the next arrives which 
is what happens in the actual situation.) Run until loss of lock and record the number of cosmic spikes 
until loss of lock. (Since loss of lock does not occur unless cosmic spikes exceeding 50 photoelectron 
counts are present, then intuitively one would figure that the probability of maintaining lock is related 

to the number of these cosmic spike events which occur.) 

b) Make 25 simulation runs varying the photon noise profile from run to run and record the 
number of cosmic ray spikes until loss of lock for each run. From tliis random sample, estimate the 
probability of maintahiing lock versus the number of cosmic splices until loss of lock. 

c) Multiply the number of cosmic splices until loss of lock by the mean time between cosmic 
spikes to determine probability versus time, 

The idea of subjecting the FGS to cosmic spikes at a higher rate than it will experience in nonnal opera- 
tion is analogous to what is known as accelerated life testing, whereby different types of components are 
tested in a far more severe environment than they will normally be expected to operate in so as to 
significantly shorten the time to failure and hence the time required for life testmg. The reliability of 
one versus another can be quickly ascertained. Sometimes, prelimmary experimentation is earned out to 
determine the relationship between the proportion of failures that can be expected under nominal con- 
ditions and under various levels of accelerated environmental conditions [3], 

The simulation conditions used in the 25 runs are similar to those used previously: 

a) PCS configuration 2 

b) Best four out of five gyros 

c) No reaction wheel induced jitter 

d) 14.5 Mv star with two axis photon noise 

e) Cosmic spikes every 2 sec in tube A, large enough to saturate the Q limiter 

f) 5400 sec simulation run times. 
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Tile simulation results in terms of the number of cosmic spikes until loss of lock are shown in 
Table 3 where the data has been ordered sudi that xj < X 2 < ... < X 25 where Xj is the number of 

cosmic spikes until loss of lock for the ith run. It should be observed that the mean and standard devia- 
tion for the random sample are 63 and 59, respectively. A histogram was constructed from the dat', to 
show the number of times loss of lock occurred (class frequency) for various intervals cf the random 
variable (the number of cosmic spikes until loss of lock). The result is shown in Figure 6 . This was then 

TABLE 3. SIMULATION RESULTS FOR THE BASELINE SYSTEM DESIGN 
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Figure 6 . Histograms based on the simulation results for the baseline system design. 
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put into tlie form of a probability density function (p.d.f.) by dividing the class frequencies by 25 (the 
total number of samples) and then dividing this by 40 (the width of each class interval). Tliis result is 
also shown in Figure 6. However, the p.d.f. derived from the random sample looks suspiciously like the 
p.d.f. of an exponential distribution which has the form 

X 

f(x) = i e ^ (1) 


where x is a specific value for the random variable X and n is the mean of the distribution of X [4] . 
If n equals the mean of the random sample or 


M = 63 


( 2 ) 


the result is shown in Figure 6. It can be seen that the p.d.f, derived from the random sample matches 
the theoretical curve quite well. Furthermore, for the exponential distribution, a, the standard deviation 
of the distribution of X, is equal to the mean of the distribution or 


a = (i 


(3) 


From Table 3, it can be seen that the standard deviation of the random sample is close to the mean of 
the random sample. This then is a second piece of information wliich infers that the number of cosmic 
spikes until loss of lock has an exponential probability distribution. A formal test of goodness of fit of 
the random data of Table 3 to an exponential distribution was performed using the Kolmogorov-Smirnov 
test.^ According to tliis test, the distribution function from a random sample is compared with that for 
a specified theoretical distribution [5]. The distribution function (d.f.) F(X), where 


F(X) = Pr(X x) 


(4) 


obtained from the random data of Table 3, was plotted along with that for the exponential distribution 
given by 


x 

F(x) = 1 - e 


(5) 


with 


M= 63 


( 6 ) 


1 . The author wishes to express liis appreciation to Mario Rheinfurth and Frank Pizzano, both of 
Marshall Space Flight Center, for their assistance in tliis area. 
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The results are shown in Figine 7. When the two curves are compared at values in the random sample, 
the absolute values for the differences in the two curves was observed. The maximum was found to be 
0.081 as seen in Figure 7. This then was compared with a critical value found in Table 4. In general, 
the Kolmogrov-Smirnov test does not apply where parameters in the theoretical distribution are estimated 
from the random sample being tested. However, the test has been worked out for the exponential dis- 
tribution wliere the mean of the distribution is unknown and must be estimated from the random sample 
at hand. Table 4 is the result of that effort [6J. A general rule of thumb is to use a level of significance 
of 0.05 (i.e., the probability of false rejection is 0.05), At that level of significance and with a sample 
size 


N = 25 , 


(7) 


it is seen that the critical value for D in the table is 


D = 0,210 


( 8 ) 


Since the value for D from Figure 7, which is 


D = 0.081 


(9) 


is less than the critical value found in Table 4, the hypothesis that the random data of Table 3 comes 
from an exponential distribution cannot be rejected. To even further argue that the number of cosmic 
spikes until loss of lock has an exponential distribution, the exponential distribution is characteristic 
of the probability disuibution for the time to failure of a system or component wliich has a constant 
failure rate [7] . Furthermore, a system or component with a constant failure rate does not wear out 
with time but is “as good as new” wliile it still functions. For the system at hand as modeled, there 
is nothing in it which wears out with time and indeed is “as good as new” wlule it maintains lock. 
Hence, by this reasoning one would expect that the time (or number of cosmic spikes) until loss of 
lock would follow an exponential distribution. 

Once the fact is accepted that the number of cosmic spikes until loss of lock has an exponential 
distribution, the real function of interest is not the p.d.f. or the d.f. but the reliability function R(X), 
which is the probability that toss of lock will occur after x number of cosmic spikes or 


R(x) = Pr (X > x) 


( 10 ) 


and is given by 


X 

R(x) = 1 - F(x) = c ^ 


( 11 ) 
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Figure 7. Distribution functions based on the simulation results for the baseline system design. 




TABLE 4, KOLMOGOROV-SMIRNOV TEST FOR THE EXPONENTIAL DISTRIBUTION WITH MEAN UNKNOWN 
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For 


63 


( 12 ) 


the reliability function, or equivalently, the probability of maintaining lock, is as shown in Figure 8 . 

However, the probability of maintaining lock versus time is preferred. This is arrived at, then, by 
multiplying the number of cosmic spikes until loss of lock in Figure 8 by the mean time between cosmic 
spikes given in Table 2. The result is shown in Figure 9. It indicates that the probability of maintaining 
lock for 10 hr (maximum observation time) under FGS control with two FGSs tracking 14.5 Mv stars is 
about 0.6 at the beginning of the mission (PAD setting = 40 photoelectron counts) and about 0.4 at the 
end of the mission (PAD setting - 80 photoelectron counts). For 24 hr under FGS control with two 
FGSs tracking 14.5 Mv stars, the probability of maintaining lock is about 0.3 at the beginning of the 
mission and about 0.1 at the end of the mission. 


V. IMPROVING THE PROBABILITY OF MAINTAINING LOCK 


The next logical question which arises is: What can be done to improve the probability of main- 
taining lock shown in Figure 9? Simulation runs were made varying two parameters in the fine lock 
loop, K 3 and both of which are fimiware parameters in the FGE, changeable in-flight from the 

ground. The various combinations of K 3 and considered are shown in Table 5. Case I corresponds 
to the nominal values for K 3 and K ^4 and hence is a reference wliich the others are compared against. 

The simulation conditions are the same as those used in the analysis in Section IV and are shown in 
Table 5, In each case, 25 simulation runs were made varymg the photon noise profile from run to run, 
the simulation run times being limited to 5400 sec in all cases. Because of the limited run times, loss 
of lock did not occur in some of the runs. However, tins can be accommodated, with the problem at 
hand, in calculating /?, the estimated value for the mean number of cosmic spikes until loss of lock. 

Since the random data is drawn from an exponential distribution, a good estimator for /a is given by the 
expression in Table 5, which accounts for the fact that some of the observations may be censored [ 8 ]. 

The variable r represents the number of runs in each case which did lose lock, (n - r) being the number 
which did not lose lock. The resulting estimates for n in each case are presented in the table. 

To see the sensitivity of to IC 3 , was plotted against IC 3 for a nominal value for as shown 
in Figure 10. It indicates that below about 0.4 for K 3 , {li is very sensitive to IC 3 and, there, a considerable 
improvement in p can be realized by a small reduction in K 3 . 

Given values for ju in each case and using the procedure followed in Section IV, the probability 
of maintaining lock versus time was determined for each case and the results are shown in Figure 11, 

The results indicate that a considerable improvement can be realized by reducing K 3 and/or K^ 4 . 
Indeed, if K 3 is set to zero, the probabihty of maintaining lock is 1 regardless of the time interval. 

The effects that reducing K 3 and K ]^4 have on other aspects of system performance need to be 
considered, but will not be addressed in this study. 
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Figure 8. Probability of maintaining lock versus the number of cosmic spikes until los of lock 

for the baseline system design. 
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Figure 9. Probability of maintaining lock versus time for the baseline design. 


TABLE 5. RESULTS OF SIMULATION VARYING K. AND K 
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Figure 11. Improving the prohahility of maintaining lock. 


VI. CONCLUSIONS 


To summaiizc liie ovci’aU conclusions or tills stuily; 
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n It was I'oniul that the time to loss of lock has an exponential probability distribution which 
is eharaetedstie of a system which has a constant failure rate. A system with a constant lailure rate 
does not wear out with time but is -“as good as new” while it still functions, Tlie FGS as modeled has 
nothing in it wliicli wears out with time and, consequently, should exhibit a constant failure rate. Hence, 
one would expect the time until loss of lock to have an exjmnential probability distribution. 


lock; 


in the baseline design, with 2 FCnSs tracking 14,5 Mv stars, tlie probability of maintaining 


ai For 10 hr is about O.b at the beginning of tlie mission (PAD setting - 40 pliotoclectron 
counts') and about 0,4 at tlie end of the mission (PAD setting « 80 pliotoelectron counts^ 

b) For 24 hr is about 0.3 at tlie beginning of the mission (PAD setting ~ 40 photoelectron 
counts) and about O.l at the end of the mission (PAD setting = 80 pliotoelectron counts). 


3) The probability of maintaining lock can be improved by reducing the values for the FGS 
fine lock loop parameters K 3 and Kj 4 , 'I'liese are firmwaw parameters in the FGE, but are changeable 

ind'liglit from tlie ground. Imleed, setting K 3 equal to zero eliminates loss of lock altogether and the 
probability of maintaining lock becomes 1 for all time. The efA’Cts that reducing IC 3 and lvj 4 have on 
other aspects of system performance need to be considered, but were not a part of this study. 
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